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ABSTRACT 


iiewecbicebave of this work was tosdecemmmnme mre expected 
deflections prior to and at the rupture load for a panel of hull plating 
of an ice~strengthened ship. An idealized mathematical model of a pin- 
nacle of ice being forced against the plating was used as the assumed 
loading mechanism. The panel of hull plating was modeled as a clamped 
rectangular plate. The idealized ice-loads were defined as a four sided 
fememid Of ice which when forced against the center of a rectangular 
Peee Crushed aft a given pressure. Thus the size of the load bearing 
surface continually increased as the applied force was increased. 


A method of analysis was established which provides a step-by- 
step solution of the large deflections of a plate and a means of deter- 
mining the minimum expected failure load. The analysis is begun at the 
Bere where the load on the plate from the pyramid of. ice is large enough 
to cause plastic collapse by a three hinge mechanism. A solution method 
for the deflection of a uniformly loaded membrane forms the basis for the 
calculations of the plate deflections. Then the energy absorbed by the 
Peoee at the point of expected failure is calculated from an equation 
that depends on the total strain at the point of failure. 


The possibility of failure of the plate at any given load and 
@eerection may be determined from the failure criterion. This criterion 
was established by forming a relation between a parameter representing 
the geometric properties and a parameter representing the material pro- 
Perties of the plate. The relation was formed from published data on 
Mewees that had been tested to failure. The calculated failure parameters 
for these plates were used to form a linear equation which correlates 
meee experimental failure loads. The criterion was then checked by 
Semeating the results given in other experiments on the large deflections 
of plates. : 


The method of analysis was then applied to plates similar to 
Poe nuli plating of ice-strengthened ships. The results showed that the 
method of analysis can be used with a reasonable degree of accuracy to 
Seculate the large deflections of a plate subject to a lateral load and 
Mempredict the load and deflection at which failure may be expected. 


Thesis Supervisor: J. Harvey Evans 
Title: Professor of Naval Architecture 
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INTRODUCTION 
Dascussion of Problem 


The solutions to two interrelated problems are sought in 
this analysis. The failure of clamped rectangular plates may be 
separated from ice-loading of plates, but the two studies are brought 
together here for the consideration of a special problem in ship's 
structures. The full analysis of hull plating subjected to the various 
loading conditions associated with navigation through areas where ice 
is present brings the two topics together. The first problem is to 
@e-wermine the ralatie tctween the hull plating and the ice that may 
be forced against the plating. The second problem is encountered as 
the force driviny the ice against the ship is increased and the hull 
plating deforms until plate failure or rupture occurs. The total pro- 
cess may be considered continuous, starting with the initial contact 
with the ice, subsequent ice crushing and plate deflection, and finally, 
extreme plastic deformation of the plate and failure. 

The need for this analysis is generated more from the design 
of merchant ships which are to operate in ice filled waters than from 
icebreaker design. JIcebreakers are generally very conservatively 
designed because they must be able to withstand ice-loadings of almost 
any magnitude during their lifetime. On the other hand, merchant ships, 
even those with icebreaker bows, have to balance their design between 
Peerit and repair costs. Their carrying capacities are greatly affected 
by the addition of thicker steel plates along the ice-belt. Furthermore, 


if there is a chance that the ship could lose part of its cargo through 


ih 





a ruptured plate, the economics of the situation increases. Even more 
mipGgistance tough, is the rising coneern over the possibility of oil 
Tolwuiien Of Lhe arctic seas. In order to prevent oil spills, it has 
been proposed that a tanker be designed with outer buffer tanks filled 
with water. The euter hull wou ldeeew@nrseninner plate than the usual 
single hull ship. The separation of the two hulls would be based on the 
maximum deflection that could be expected from the outer plating when it 


failed under ice-loads. 
Tdealization of Ice-Loads 


ane ice-loads to be considered will be idealized with regards 

to their geometry and point of application. The present design practices 
are based on static loading of the hull plates. The model loads used 
here will be considered continuous but applied in incremental static steps. 
Icebreaker and ice-belt plate scantlings are usually obtained by using a 
uniformly i oGatouincd pressure Equal tOurbe 1¢Ge criushine pressire. Since 
the actual ice crushing pressure is not constant and dependent upon 
geographical location, age, salinity and many other factors, it was 
desired to try to cover a wide range of possible loading mechanisms. 

The driving force behind the plate/ice interaction must also be considered. 
It was for these reasons that a pyramid shaped pinnacle of ice was chosen 


as the idealized ice-loading mechanism. 


Type of Analysis 


Plate failure analysis lies between the theory of plasticity 


and a handful of experimental tests that have carried applied loads up 


ee 





Eour dtecmmupeune. Theoretically it would be desirable to use the theory 
SBreplacticity Lo analyze the plate deflection and failure. This was 
Pemeicgered, but after reading other works in the field and considering 
Ehe accuracy of the possible results, it was decided to use a more 
empirical approach. The final analysis will be a combination of plastic 
analysis and engineering judgment. This combination of approaches to 

the problem allows a study of the plate/ice interaction within a limited 


mimes period with sufficient accuracy of the results. 
Load/Deflection Relations 


Most of the By pe cieneall work on rectangular plates with large 
deformation has been done under uniformly distributed loads. To obtain 
any data on concentrated loads, it was necessary to go back to a report 
femerperiments conducted in 1937. The plate deflections and stresses 
Gemeekept in the elastic ree in these tests. The results did provide 
a good starting point for the shapes of the plates during plastic defor- 
mation under concentrated and partial loads. Some of the results and 
observations on circular plates under concentrated loads were borrowed 
for the analysis. 

It was necessary to consider the concentrated or partial loads 
in order to present the complete model of plate/ice interaction. The 
crushing process of ice in contact with the ship's hull will begin as 
eeeoncentrated load. The ultimate strength of sea ice in compression 
is usually low enough in relation to the design pressures of hull plating 
Pu@emechat the ice-loads are uniformly distributed before permanent defor- 


mation of the plate begins. However, since the ice crushing process will 


r 





require some partial loading, large plate deformations may occur under 
Pravlenr.spuL more concentrated loads. Eherefore, partial loading of the 
plate will be analyzed in order to consider all possible extremes for 


the failure analysis of the plate, 
Failure Analysis 


Experimental results of rectangular plate failures are usually 
an added column of data to the experimenter's primary reason for testing. 
The beginning of investigations into large plastic deformation and failure 
imemerecen attributed to Bach. His experiments conducted in 1908 formed 
the basis for Hovgaard's 20% permanent deflection rule.! Messe Cesics. 
like most of the data obtained for rectangular plate failure, were for 
f@eormly distributed loads. Often the comment that the plate failed at 


(14) in 


a given pressure is the only amount of data given. Greenspon, 


: E2,E3 
Sweuuating experiments by J. W. Day, a) attempted to correlate the 


—-_ = 


Failure data so that plate rupture under a uniformly distributed load 
could be predicted. His recommendation was to establish a limiting 
maximum deflection to span ratio. Using a statistical approach, it could 
be stated that a uniform lateral load giving a deflection to span ratio 
Geeeeeer than 0.10 could be considered approaching the plate failure 
load. The only other author who discussed plate eee any extent 
was Loeser. — 

Most of the experiments as reported by the authors listed in 
tne Dibliography under Plate Deflection-Experimental were helpful in 


establishing a failure criterion. It was decided early in the analysis 


that an energy relation would be used. This enabled correlation of the 


1 See reference T5, pp. 120-125. 
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failure data that was available. Not all of the data was useful, since 
all material properties must be known in order to properly correlate the 
experiments. Since Day and Loeser included all material properties in 
muerte ports, the results of their experiments were relied on most 


heavily. 
Outline of PiEOced tine) 


The procedure of investigation that follows will continue in 
a similar format as used in the introduction. First, the model of the 
ice-loads will be formed. Then, the plate deflection/load relations 
will be established. ne each increment of loading or deformation of the 
plate, the amount of energy absorbed by the plate material should be 
calculated. The relations for this segment of the analysis will depend 
on the load/deflection calculations and will serve to relate them to the 
Pwoiee criterion, The Final Ste peli Ehemanadiysic will bette estab- 


Meeument Of the failure criterion. This will complete the model of the 


ice/plate interaction. 


Ordering of the Bibliography 


It has been previously mentioned that this analysis falls 
between the basic theory of plasticity and empirical relations. Many 
of the assumptions and choices of methods are based on the experiments 
and observations of other authors. Each phase of the analysis can be 
Pemeeormed by other techniques. However, within the time limits of the 
investigation, all cases were carefully studied so as to make use of the 


best method available. For these reasons, it was felt that the list of 


HES) 





references should reflect the degree of influence of each work upon the 
final analysis. Therefore, the references within each subject area are 


listed in approximate descending order of influence or use. 
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PROCEDURE FOR CALCULATING ICE/PLATE INTERACTION 
Initial Assumptions for Ice-Loads 


A four sided pyramid of ice as shown in Figure 1 will be the 
loading mechanism. The manner of loading will be such that it may be 
@omsrdered to be static at all times. 

The included angle of the pyramid frustum will not effect the 
load bearing area ae ice on the plate. Also it is small enough to 
eamee only the crushing portion of the ice to make contact, That is, 
maemo oped sides of the cone will not come in contact with the plate at 
any time. 

The load bearing surface will conform to the curved shape of 
Biemplate as if the plate in the region of the load is subjected to a 
vniformly distributed load. Thus, load edge stress concentrations will 
be neglected. Also, at no time will the area of load bearing earn be 
tess than twice the thickness. Therefore, shear forces will be neglected. 

The load bearing area will be a square of side equal to C. 

If C equals the width of the plate, the load bearing area will be set 
equal to the plate area. Also, at this Pint, the loading pressure will 
fe eomieted from a2 constant ice crushing pressure to a pressure equal to 


the total load divided by the plate area. 


Initial Assumptions for Rectangular Plates 


The plate material will be isotropic and homogeneous. All 
four sides, of the plate will have a clamped edge condition for support 


as shown in Figure 2 and only in-plane stresses are to be considered. 
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Figure 1. Idealized Ice-Load 
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Only deflections after the plastic collapse load is reached will be 
considered. A rigid, perfectly plastic material is assumed for the 
collapse mechanism and a biaxial strain condition will be assumed at 


the mid-point of the plate. 
Model of Ice-Loads 


The idealization of an ice-load will be a four sided pyramid 
of ice. The actual applied load will be in the form of a force acting 
Gomene base of the pyramid. The force will be transmitted to the rec- 
tangular plate in the form of a uniformly distributed load acting over 
an area centered at the mid-point of the plate. The value of the load- 
fmemeprescstre will be the ice crushing pressure except during initial 
loading and after the bearing surface is equal to the plate surface 
@rea. These exceptions will be discussed later. 

Aomene Grivime force increases, the bearing surface of the 
pyramid will increase as the ice crushes. Ths a dt Cerne wT Cemerusil Tae 
Peeesure is considered constant throughout the loading process, the 
Deereime surface area will increase directly with the driving force. 

Since the area, e Semiiy einen is equal Co the driving foree, 
F, divided by the ice crushing pressure, Ee Siewlemebhmor tencldemcat 
the bearing area will be: 


ga (1) 


fiterseetelation will hold after the pyramid first starts to crush until 
the plate is completely covered by the bearing surface. The loading 


meer to the initiatiom of ice crushing will not be considered, since 
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the plate deflection will be in the elastic range. 

Once the bearing surface equals the area of the plate, the 
applied load, F, will remain constant unless Fi is allowed to increase. 
Thus it will be assumed that ice crushing ceases and P will increase 
with any increase in F as; 

Pp = F (2) 


Plastic Collapse Load 


The collapse mechanism for the plastic anaeene is the general 
memuesnaped deflection surface. For the transverse cross section a 
three hinge collapse mechanism is used. For any loading, it was assumed 
that the center plastic hinge forms either simultaneously with the two 
edge hinges or at some increased load beyond the formation of the edge 
hinges. = 

Since the analysis is for large deflections of plates the 
deilection at the plastic collapse load will be the first to be considered. 
feeeertunately, there is no easy solution for the deflection at this load. 
The problem is further complicated when partial loads are considered. 
Peee displacement relations from the theory of plasticity and a 
membrane solution were considered to be the easiest methods for finding 
Meemeentections. A similar degree of accuracy asi be obtained from 
either method. The membrane solution has been shown to provide good 
representation of the plate deformation for loads well beyond the collapse 
load. Therefore, it was decided to use the membrane solution throughout 


the analysis. 


Zl 





As the pyramid of ice is forced against the plate, the pressure 
Peeinemover the initial bearing surface will rise until, the ice crushing 
pressure is reached. During this time, the plate will behave elastically. 
From the time pe is first reached until plastic collapse of the plate, 
the area of the bearing surface will increase. There will be associated 
with the collapse load Fy? a constant Pe and an area is inrordeta ic 
use the membrane solution to find the deflection at collapse, the value 
OL E.. must be known. Actually it is easier to use the non-dimensional 
form C/A. 

The moment M. at the center of the plate was found using an 
elastic solution for clamped plates under partial loads. M.. = e’ P 
where P is the total load and the load coefficient 6’ is given in 
Figure Al for values of B/A and C/A. Then by setting M. = Mo = Ono /4, 
Meemeaela moment forsa rigid, perfectly plastic material, the bearing 
surface area can be related to the yield stress and plate thickness. 
Thus, when P = pc? and both sides of the equation are divided by the 


width of the plate, A, the following non-dimensional equation is 


eo *L (a) (a) 7 


From equation (3) the bearing surface area at the collapse 


formed: 


load may be determined for any plate. Various solutions to the equation 


are given in Figures A3-A5 for given values of B/A and A/H. 


Z This is equation (A7) of Appendix A. 
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Plate Deflection Calculations 


Once the size is determined for the load bearing surface at 
the collapse load, it is necessary to calculate the deflection of the 
plate at this load. As discussed previously, a membrane solution will 


Memused for the deflection calculations. The membrane formala,” 


7 4qA7k (4) 


in general only holds for a uniformly distributed pressure over the 
entire plate. This meant that some geometric relationships were re- 
quired in order to use the membrane formula with partially loaded 
mieies. It was assumed that the plate in the region of the partial 
Bitesc behaves as a membrane after plastic collapse was reached. Then 
the membrane or in-plane forces along the edge of the loaded area are 
Carried through to the edge of the plate. This implies that the slope 
of the deflection profile at the edge of the loaded area is the ee 
as the slope at the edge of the plate. Thus the conical shape is 
represented for the concentrated loads and the full membrane shape is 
Meeaewhen the load is uniformly distributed over the plate. It is 
Beweved that this provides a good representation of the true deflection: 
The deflections for the concentrated loads may be less than they should 
beeebut the accuracy should improve as the load bearing area increases. 
The total deflection calculations consist of three basic 


Peepers rirst, the deflection Z for the loaded portion of the plate is 


found from the membrane formula. For the width of the applied load C 


See Appendix B. 


See figure 3. 
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and the pressure oe Dicureeomemevycuthe relation of Z to the total 
@eelection. The second step is to determine the- slope of the plate 
ats Bie pace OL the foaded area. fhe PP aciue of the deflection 
profile provides the slope at the point x = (Af C)/2. Finally, the 
total deflection at the mid-point of the plate is found by combining 
the membrane solution with the assumed slope at the plate edge. A 
eemple geometric relation will determine the deflection at the edge 
of the loaded area, and the membrane solution may be added to it to 
eive the total deflection. The central deflection W is given by: 


W=Z+1 (A-C) tany , : (5) 
5 


where the membrane formula 


Z=4 c He (6) 
™ 


Po the detlection of the Loaded portion, 


BueEan wee G2 = tee eka atx =A™7C, Cp 


dx Zo 2 





The constants k and ka Agewmthes factions that Bepend sere... Lhe wages 
for these constants are included in Appendix B for given plate aspect 
ratios. 

Prior to discussion of the procedure for load/deflection 
analysis, the membrane solution for plate deflection must be understood. 
mmeeeeecual formulation of the deflection relation is given in Appendix 
Meeiaere are two dependent variables in the relation. The membrane 
Perce , S, which acts along the edge of the loaded area depends on the 


total deflection of the membrane, W, which in turn, depends on S. The 





2 


See Figure Bl. 





usual procedure for obtaining the value of S is to assume a membrane 
deflection and shape. Then a strain relation is used to obtain the 
principal strain in the plane of the plate for the assumed deflection. 
The usual assumed deflection shapes are cylindrical or parabolic. 
invordem to model themice Voedincsconeuemon.) the equation for 
the shape of a partially loaded, simply supported plate was used. A 
biaxial strain relation which depends on the B/A ratio, the size of the 
bearing surface given by the C/A ratio, as well as the maximum deflection 
of the plate was established for the mid-point of the meen This made 
it possible to obtain the membrane stress from the stress-strain curve 
for the plate material, Ths, Lor a given deflection, the membrane 
fo inee Mit length is equal to the membrane stress times the plate 
thickness. 

If the membrane formula (equation 4) is used alone, the applied 
load, q, could be found that would produce an assumed deflection. Since, 
for the model of ice-loads, the loading pressure is assumed constant, 
pumeeerative procedure must be used to obtain the proper correlation 
between Z and S for a given constant load lee 

For the calculations involving the deflection of a partially 
femded plate, the total deflection, W, is used instead of Z to determine 
Ene value of the membrane force, S. This appears to provide a more 
realistic situation, because the strain relation as derived in Appendix 


Deadepends on both the plate shape and deflection. 


See Appendix D. 
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Load/Deflection Relations 


Tne procedure for sopra imine ene FeTee rion Wipe Obmeaen we ive 
of the total load, F, will be iterative. The beginning will be based 
on the plastic collapse deflection. After the value of (C/A) | is found 
from the plastic collapse relation, the membrane solution will be used 
Bomrind the deflection at the plastic collapse load, E Alle Geabetete 
approximation for the deflection is found by assuming that the membrane 
stress is equal to the yield stress. Then for the value found for Co 
from (C/A) | and the ice crushing pressure, P,, the first values of Z, 
dzj/dx and W (equations 6, 7 and 5) can be calculated. This ie the 
starting point OE ythe iteration for the proper value of Wo for the given 
load FO. The successive solutions of W will follow the same format as 
outlined below for the general solution. 

For any given total load By there will be a corresponding 
bearing surface area represented by the length of one side C_. The 
Peeveem is to find the deflection W, that Corresponds G0 es and Co 

Let FO = Fo [Toecy where the supccripe n designates a certam 
value of static load, and AF is the incremental value of the total 
toad. The subscript o refers to the previous value of applied load. 
For example, an Vi intel aiverenresemt thescollapsesload, Then 


Ce ee erie (8) 


will provide the necessary load bearing area to determine the value 
or Ww. . 
i 


The previous value for membrane force will be used to obtain 


vay 





a membrane deflection for C. ily ee ae Issa 


By th Goa te Ps | (9) 


mage £oLlows that 


(GZ = Pe Gn _ kd (10) 
dx i Zs 
and 
W. = 4,4 | (4-C,) az Gil.) 
1 7 — — , 
2 dx 1 


The new values We and Cs are used to enter the biaxial strain rela- 


tions and obtain a new value for the membrane force Si. Theses vse p 


memeerreturn to equations (9), (10), and (11) to calculate Zea: 


(dzZ/dx) and Weaye The subscript i+l is used to distinguish the new 


oe” 
Memes Of deflection from the membrane force. This notation indicates 
that the deflection calculation is made with the membrane force found 
meomeene previous deflection, W,. After each new value of W. is found, 
its corresponding membrane force S. is determined from the strain 
relations and compared with the previous value of S. Thus, the calcu- 
lation will be concluded when the two values of the membrane force 
compare favorably. 

iremapoverprecedire 15 repeated for each new value of ff. 
When the load is reached which produces a load bearing area equal to 
Siemarea Of the plate, that is when C = A, then the loading pressure 
Ba feel cease to be constant and increase with the total load. For 
met se quent increases in load C/A = 1.0 and P, = F,/(A x B). This 
assumption of continuous loading on the plate is just an extension of 
the model of the ice loading beyond the point where the pyramid of ice 


bes been fully crushed. 
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Energy Absorbed Relations 


Once the load/deflection relation is found, Ene energy 
absorbed by the plate should be determined at each value of F in order 
to evaluate the possibility of plate failure. The energy of deflection 
fomeimbe Che greatest at the edge of the plate except for cases of very 
small load bearing surfaces. For the plate thicknesses usually con- 
sidered for ships, the load bearing areas will be large enough to rule 
Out the shear effects of the comeentrated Loads, 

Sincemtie hate Faibure Can be expected at ane mid-point of 
the long side of the plate, the energy absorbed at this point should be 
Sreater than that at any other point. There will be two sources of 
strain energy at the edge of the plate. The initial source will be from 
the bending moment. Then as the membrane forces increase, the portion 
of the total strain due to the bending moment will decrease. The com- 
bination of the strain due to bending and stretching will give the 
total strain at the expected point of failure. 

Using the assumption that plane cross-sections remain plane 
during the bending and stretching, the strain distribution across the 
plate thickness will be as shown in Figure 4. The rotation of the 
cross-section is due to the plate bending, and the location of the 
neutral axis away from the mid-plane is earned bymene ‘Siz erotic 
membrane stress. The influence of the membrane stress is equivalent to 
lowering the neutral axis from the geometric center to a new point 2}. 


Tots shift in axis of rotation of the cross~section depends on the 
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PLATE MID -PLANE 


Figure 4. Strain Distribution for Combined Bending and 
Membrane Forces 
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ratio of the membrane stress to the equivalent yield er cee 


2 0, 
where 0, = 0, See OF home =O. 5 (13) 


V1i- V2 p2 


ieeesrelation will ho:d true for values of OF greater than OO since 
femenat point the neutral axis will lie “outside” of the material and 
all strains across the thickness will be in tension. 

ine Strain Of the outermost Tiber al Ehe mid-point of the 


long side will be represented by Cc .- From Fipure 4, 


{} 


ei. +e tan 86 (14) 


Cz 


Since E . = 2, tan 8 = : O. tan 6 


Op 


ll 


C. H tan 6 Oa (15) 
Z Op a 


Most experiments have shown that th: longitudinal strain, 

ae at the mid-point of the long side of the plate is very small 
Compared to the transverse strain at that point, Tae on was taken 
as the principal or transverse strain at the point of expected failure. 
Again using a biaxial strain relation as given in Appendix C, the 


uniaxial strain is, 


oF . ay (16) 


where V is equal to Poisson's Ratio. 


Young in reference E5 recommended the use ot JQ) in the equation 
for the shift of the neutral axis due to the ildition of membrane 
Bene sses. 
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[MimondcratO Lindethe stpesseat the expected point of £Eailure, 
the uniaxial strain was used in the constituative relation for the 
Peete material. The total stress will be designated o and may be 
found for the strain €, from the stress-strain diagram of the material 
in simple tension. 

The energy absorbed per unit volume at the mid-point of the 
long side will then be represented by: 


ie On fincaos | Gey) 


eta 
This value may be found at any point in the loading sequence if it is 
assumed that 6, the angle of rotation of the cross-section, is equal 
meme, the angle of the slope of the deflection profile at the edge 
of the plate as shown in Figure 3. This is a reasonable, as well as 
a convenient assumption, since the slope at the edge is determined 
Gmeine the deflection calculations. Thus, for each value of F and its 
Peeaondine W, 


tan 6 = dz (18) 
dx 


neva , may bem Gewritteml as 
esate (2 }( 2] (19) 
gia) Waa 
Also, the stress-strain relations for the plate material that are 
required for the deflection relations may be used to find O |: There- 
fore, Lor any loading and corresponding deflection, the energy absorbed 
Hyetne plate at the mid-point of the long side may be determined. A 


comparison of this energy to the restraint energy of the material will 


fee to a failure criterion. 


eZ 





Failure Criterion 


When a plate is drawn in biaxial tension, the strain rate 
becomes unstable when the principal stress reaches the ultimate strength 
Bieene material. The plate begins to decrease in thickness rapidly like 
a tension specimen begins to neck at the ultimate eee For this 
reason it was decided to use the ultimate stress of the plate material 
fo establish the failure criterion. The ultimate stress will be defined 
as the stress at which the slope of an engineering stress-strain diagram 
memedual zero, The area under the stress-strain curve for all strains 
feeewtnan the strain at the ultimate stress will represent the stability 
of the plate material. 


A restraint energy relation will be defined as 


R= € O . Ey (20) 


a aT 
where O , is the ultimate stress and 7 is the ultimate strain. 
Both values are defined at the point of zero slope of the stress-strain 
euinye . 

Hheweetmm restraint enerey is used Since R represents the 
amount of energy that might be absorbed by the material before insta- 
Deeiity begins. A more aecirate measure of the ability of the material 
to absorb energy would be the area under the stress-strain curve. Since 
Bueomonea is usually difficult to obtain readily, the R relation is used. 
It provides a reasonable representation eee auea undememost stre@ss- 
strain curves. Furthermore, it is more desirable to obtain a relation 
that will give a consistent measure of the restraint energy rather than 


Gmemoccurate one. 
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ihewiwo paramerers that bescn correlate the minimum amount 
meenarlure data that is available are the deflection to- span ratio, 

W/A, and the absorbed energy to restraint energy relation U/R. These 
non-dimensional parameters represent virtually all of the factors 
involved in the plate problem. Several other parameters were considered, 
but there did not appear to be any relation between them and the failure 
data. The W/A ratio helps to bring together all the geometric proper- 
ties of the plate. Tees eal properties are well represented by the 
ine rat Lo. 

Figure 5 is a plot of the two non-dimensional parameters for 
most of the rectangular plate failure data that is available. The curves 
are drawn through data points obtained by the calculation procedure 
m@eeerrped in the previous sections. The failure points for each plate 
are designated in two positions. ine calculations for the det leer. 
Gemenie plate did not correspond directly with the experimental jae 
femme references E2 and E3. In general, the calculated deflections 
were slightly Weve etnan the experimeneal deflections for a given toral 
load. The possible reasons for this discrepancy will be discussed later. 
ieee telt that both failure deflection and failure load should be 
indicated in Figure 5. Thus, the first set of failure points were found 
meres the experimental deflection at the mid-point of the plate to 
determine the W/A ratio at failure. The second set of points were ob-~ 
tained by using the failure load from the experimental data. For each 
rupture pressure, the total loading force F was calculated and from 
Figure 8, the corresponding U was found. The U/R values for the failure 
loads are shown on Figure 5 as the lower set of failure points on the 


Curves. 
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The lower set of failure points was chosen to represent the 
failure data. Although the lower set provides the more conservative 
choice for a failure criterion, the major reason for using this set is 
that the failure data corresponds better to the deflection and absorbed 
energy calculations described in the previous sections. 

Ae tine can be drawn on Pigumeso connecting amos Of /themratl— 
ure points. With the exception of the plate designated 1B this line 
represents a minimum relation of W/A to U/R for plate failure. If it 
is moved slightly to the left, the failure point for plate 1B will also 
lie above the line. This new position of the limit line will be used 
as the failure Pte LTIGN, 

If §& = W/A and @ = U/R, then failure may be expected beyond 
the line defined by 

Ome one. 12.0 (21) 

Greenspon in reference T4 used essentially the same failure 
Bea and stated that failure should be expected for loads beyond the 
pressure that produces W/A = 0.10. This limit as shown in Figure 5 is 
quite conservative for 6 <0.9, but for plates with 6 < 0.10 and > 0.9 
the limit may be too large. 

The energy relation was chosen for this analysis in order to 
Meeovede a parameter that will accurately represent all ae loading 
conditions. Unfortunately, the available failure data is for loads 
uniformly distributed over the entire Sees inporder born uope muy 
relate the failure for partial fonds and uniform loads, equivalent 
absorbed energy relations are used. Thus the failure criterion should 


hold for all but concentrated loads with small bearing surfaces. 


o> 
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Summary of Method Or Analysis 


The method of calculation of the deflection and the absorbed 
energy for any plate has been incorporated into a short computer program. 
wieewaescription of the program itself is in Appendix E. The only dif- 
ference between the procedure summarized below and the computer program 
is the method of determining a stress for a given strain. A generalized 
equation for the plastic portion of the stress-strain curve is used in 
the program. With the proper choice of variables the generalized stress- 
Seen equation can represent the actual stress-strain curve of a given 
material with reasonable accuracy. There are several variations to the 
general equation that may be used. The relation described in Appendix 
Beappeared to be one of the easiest to use. 

ihewecOnpUuceieDLovided Chesspecd and =aceliLacyabtoucal ry) Ole 
Meemealculations for many plates. The hand calculations are not diffi- 
eult, but they are quite tedious. With the use of the tables and curves 
presented in the Appendixes, the time for calculation can be shortened 
considerably. The step-by~step calculations for a plate of known geo- 
metric and material properties and a given ice~crushing pressure are; 

1. For the values of B/A, A/H, and O76/P. enter the curves 
given in Figures A3, A4 or A5 to obtain the value of C/A at the collapse 


load. Then calculate Ne & C7 Poe 


402k Pe , and dZ = PrCkg, 
hee dx 2s 


Pemecatculate Z 


and W = Z+ 1 (A-C) _dZ where k and kg may be 
dx 


bo 


obtained from Figure Bl for B/A. Also the membrane force for the first 


oy 





approximation will be assumed to equal the yield stress times the 
thickness. 

3. The quantities W/A, B/A, and C/A are required to find 
the value of the average transverse strain at the mid-point of the plate. 
Tables of values of the average transverse strain and representative 
curves are given in Appendix D. 

4. The new value of membrane stress will be found from 
Bmesstress-strain curve for the plate material for the average trans- 
Pee strain found in step number 3. Then S; = Oe sae 

». Steps 2, 3, and 4 are now repeated as described on page 
28 until the new sities of S; is within an accepted difference from the 
previous value. (An absolute difference of 25 psi is used as a tolerance 
in the computer program). 


6. For the final values of Og and dZ/dx from step 5, cal- 
culate €, = H dz Ce Fea 


7. Again enter the stress-strain curve for the plate material 
to find on which corresponds to €,. 


Be icalcutate the eneresy absorbed at the point of expected 


failure. U 


Er, O+ 
9. Use the value for €,, as defined on page 33, and calcu- 
foen= ©, O.,. 


10. Using §& = W/A and @ = U/R, solve the fan ine Cr ernon 


pe toe 0 = 2275, 12 FE*= OO, then failure can be 


it 


equation FC 


expected for any additional loading. 
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Mein he = ©, Sinerease the total, toad F by the) load incre- 
ment A F which may be any predetermined amount. Obtain a new value 
for C = VF/P. 

[eee poate Ste pS MEhGeucimietO find tnemacs Lectionsand 
absorbed energy values for the new loading condition. 

It is not necessary to evaluate the possibility of failure 
memeach load increment. After one or two times through the procedure 
the approximate point of failure may be apparent and the solution 


method can be shortened. 
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RESULTS OF APPLICATIONS OF THE METHOD OF ANALYSIS 
Plastic Collapse Load and Deflection 


The method for determining the total load for plastic collapse 
depends on a given ice crushing pressure, The area over which this 
Peecsure must act to cause collapse is determined as outlined on page 22. 
If this method is reversed so that for a given C/A and A/H a certain 

oP. is found, then the required pressure for collapse will be known. 
Witeme this pressure is compared to elastic experimental data for plates 
of the same size and material properties, the collapse pressure should 
be considerably greater than the elastic loading pressure. It is very 
difficult to compare deflections, although the deflection should be 
increased by an amount similar to the loading. 

The plastic collapse load and deflection were calculated for 


— 


Eeitees Of the same size and material properties as those used by Sturm 


(E4) 


and Moore The size of the bearing surface was fixed at 2" x 2" 
and 12" x 12". Then the plastic collapsé load was found from Figures 
A3-A5. In this case P,, which has been called the ice crushing pres- 
Seeeeets the collapse load. Finally, using these values for C and P, 
Eremeotlapse deflection, Wo, was determined by steps 2-5 of the previ- 
ous section. The results of the calculations for three of the plates 
are shown in Table 1 along with experimental data from reference E4, 

A second comparison of the plastic collapse loads and deflec- 
tions was made with the results from the large deflection theory of 


plates. This provided a means of judging the accuracy of the membrane 


solution method for calculating the plastic collapse deflections for 
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uniformly loaded plates. 

Reference E8 compares many of the solutions to the large 
deflections of plates to experimental results. Most of the relations 
are given in the form of curves for the non-dimensional deflection, 
‘stress, and pressure ratios. The deflections of three plates were 
calculated from the corresponding pressure ratio. The plates were of 
the same size and material as those given in Table 1. The collapse 
pressure obtained from Figures A3-A5 for C/A = 1.0 was used in the 
Seessure ratio, 


eRe poAt (22) 
*  16EH4 


For example, Figure A3 shows that QO ./P, = 215 for an aluminum 

(E = e710 psi) plate 1"x48"x48" with oe = 38800 psi and C/A = 1.0. 
Thus Be Po = 180.5 psi. Using the value in the pressure ratio, a 
deflection ratio W/H = 0.9 was obtained from Figure 24 of reference.E8. 
The same procedure was followed for the other two plates. The results 
from the large deflection theory and the membrane solution for collapse 


deflections and loads Pe shown in Table 2. 
The Membrane Solution for Plate Deflections 


It was relatively easy to determine the accuracy of the mem- 


brane solutions for the plate deflection. The plate dimensions and 


E2 ,E3 
: ( ) 


material properties of the test plates used by Da and Loeser 


1 , 
=) were put into the computer program to determine the load/deflection 


Peremtenships. The results of the calculations as well as the input 


data are shown in Tables El to Ell of Appendix E. It should be noted 
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Ghat the input variable for ice crushing pressure for these calculations 
is only the loading pressure that is used for the first.calculation. 
The pressure used for any additional loading is equal to the applied 
force divided by the plate area, [It was necessary to follow this pro- 
Peaqure because all of the test data is for uniformly distributed loads. 
Actually the solution method is the same as if an ice-load is used and 
C/A = 1.0. In this way the computer program may be used for any type 
of lateral load except highly concentrated loads. 

Also, the first calculation for deflection given in Tables 
El to Ell does not represent the collapse deflection, The starting 
point used for these calculations was chosen for convenience in com- 
Peamene the results to the respective experimental data. For example, 
the rupture pressure given in Day's experiments was divided by five and 
the result was used as the input for ice crushing pressure in the com- 
P@eereprogram. Then the increment of applied force AF was set Sail 
to the ice crushing pressure times the plate area. 

The results of the calculations are shown-in Figures 6 and 7. 
The experimental results of Loeser and Day are also plotted in each 


figure. 


Energy Absorbed Relation 


MErwas ©f interest to plot curves’ of total load F to the 
‘energy absorbed U for three sets of data. The F versus U relations 
for the calculations on Day's test plates as given in Tables El to E7 
are shown in Figure 8. Similar curves for calculations on Loeser's 


test plates (Tables E8 to Ell) are plotted in Figure 9. The results 


Ad 





of calculations made on arbitrary plates under idealized ice-loads 
are given in Tables El2 to E19, and the F versus U curves are shown 
meee seure 10, 

The calculations for Loeser's and Day's test plates have 
‘been previously described. The results shown in Tables El2 to E19 
are for four plates with the same B/A ratio, but of varying thickness. 
Also, two sets of material properties were used for each plate. The 
calculations follow the format as described in the method of analysis 
pmamappendix E. 

A particular case is encountered in the calculations for 
the arbitrary plates eae was bypassed for the calculations of Loeser's 
and Day's test plates. The assumption that the ice crushing pressure 
changes when C/A = 1.0 has been discussed previously. The analysis of 
the arbitrary plates required this change. This on be seen by traci. 
down the C/A column in Tables E12 to E19, For all but the cree : 
Rees @ sostarts at the value for (Ci Ae and increases until C = A. 
For square plates the ice crushing pressure will change to the total 
load divided by the plate area. But for plates where B/A > 1.0, this 
change can cause a decrease in the crushing pressure rather than the 
desired increase. As a consequence the deflection will be smaller than 
the previous deflection. In order to help eliminate this unnatural 
discontinuity, it was assumed that the loading surface continued to 
grow after C/A = 1.0. That is, the loading pressure was kept constant 
Paci) not just the width of the plate but the entire plate was covered. 
Then for any further increases in the load, the pressure increased as 


the force divided by the area of the plate. This scheme would be the 
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natural procedure to follow if the calculations were being done by 


hand but it required a few extra statements in the computer program. 


Application OF Failure Criterion 


The relation between § = W/A and @ = U/R has been plotted 
in Figure 5 for the results shown in Tables El to E7. The failure 
criterion was established from these curves and the failure points 


Spee | 
y ( ee The establishment of the & versus ¢ curves is 


given by Da 
discussed on pages 33-35. 

For the W and U values given in Tables E8 to Ell, the corres- 
ponding §&'s and @'s were calculated. The relation between the two 
failure parameters is shown in Figure ll. These curves represent the 
experimental work of Loeser aD) 

Similarly, Figure 12 is a plot of & and @ for the calcula- 
tions made on the arbitrary plates given in Tables El2 to E19. The 
Peerary plates by their sizes and properties represent typical hull 
plating on ice~strengthened ships. Thus, the curves of Figure 12 


give a good indication of the possibilities of failure of these plates 


under the idealized ice-loads. 
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PLATE NOS FROM REFS. E2 AND E3 
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Figure 8. Load versus Absorbed Energy for Day's 
(E2, E3) Experimental Test Plates 


49 








PLATE NOS, FROM REF. £1 
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Figure 9. Load versus Absorbed Energy for Loeser's 
(El) Experimental Test Plate 
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Figure 10, Load versus Absorbed Energy for Arbitrary 
Plates of Tables El2 to E19 
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DISCUSSION OF RESULTS 
Load/Deflection Relations 


Figures 5 to ll and Tables 1 and 2 are presented so that a 
comparison of the various portions of the method of analysis may be 
made with associated experimental data. The degree of accuracy of each 
step of the analysis is shown in the tables and figures. 

Although the data presented in reference E4 was published in 
1937, it proved to be the most useful experimental data on rectangular 
plates under partial loads. The comparison between Sturm and Moore's 
experiments and the calculated values for collapse deflection and load 
Given in Table 1 is Veimywcrudenm= lhe experimental data is for the 
elastic range only, and the calculated values assume that the plastic 
collapse load has been reached. However, the experimental data does 
provide at least a lower bound to the calculated values. 

There are two general points that should be recognized from 
Table 1. The relation between the deflections and the B/A ratio that 
is indicated by the experimental data holds for the calculated values. 
freis, the experimental deflections increase as B/A increases and 
the same holds true for the calculated values. A similar comparison 
between deflections and plate thicknesses shows that the experimental 
derlections increase as the plate thicknesses decrease. However, the 
calculated deflections do not follow the same trend as the experimental 
results. It seems that the thinner plate should have the greatest 
deflection as it does in the experimental data. Although only relative 


comparisons have been made, the failure of the membrane solutions for 


54 





the collapse deflections to compare favorably with the experimental 
results should be considered a limitation of the method of analysis. 

A comparison of the plastic collapse deflection obtained from 
the membrane solution with the deflections calculated from large deflec-~- 
tion theory may be made from Table 2. In general, the membrane solution 
deflections are lower than those obtained with the same pressure from 
mamee Ceflection theory. For this case of uniformly distributed loads, 
both sets of calculations show that the deflection increases as the B/A 
decreases. Also, the deflection for the large deflection theory remains 
constant for the three plates considered. Again the relation between 
deflections and plate thicknesses does not compare eavorablye The 
membrane Petotion maintains approximately the same W/H ratio at collapse 
meget thicknesses. I[f the results from the large deflection theory 
are used to define the plastic collapse condition, then the limitation 
in the method of analysis as discussed for Table 1 is again encountered, 

The accuracy of the membrane solution may be easily seen in 
feeres 6 and 7, The calculated deaiecaer are either completely above 
the experimental deflections or completely below. In all cases, if the 
plate has strength properties (see Tables El to 11) greater than mild 
steel ( Os = 38000 psi) the deflections lie below the experimental 
data. For the higher strength plates, the deflections are all below 
the experimental results. The one exception is Day's test plate 2A in 
Figure 7. It has a yield strength slightly above the mild steel plates 
Me che curves for the calculated and experimental deflections coincide. 

Since the shapes of the calculated and experimental deflection 


curves agree very well, the only major discrepancy is the difference in 
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the deflections as noted above. The reason for the differences is the 
solution method for obtaining the membrane force for the membrane for- 
mula. This same reason may be cited as the cause for giving lower 
values of deflection for high strength materials and higher deflections 
for lower strength materials. The generalized stress-strain curve used 
in the computer program usually gives values of stress that are lower 
than the actual value for the same strain. This will in turn produce 
larger membrane deflections since the membrane force is S = O .H. 
Furthermore, for two plates of approximately the same thickness, the 
membrane stress will be greater for a plate made of high strength metal 
than one with a lower yield strength. Consequently, the membrane de- 
flections for the high strength material will be less than those of 
materials with a lower yield stress. 

For the general discussion of the application of the method 
of analysis, Figures 13 and 14 illustrate what the Wenerice eae 
relationship is for hull plating under the idealized ice-loads. The 
curves and the respective results of the computer calculations given 
in Tables E12 to E19 show that for small load bearing surfaces the 
collapse deflections appear to be excessive. the collapse deflections 
for the large bearing surfaces are usually smaller than what might be 
expected. For the uniformly loaded plates, the collapse deflection is 
not necessary for the load deflection analysis. Thus, for plates that 
require a collapse load greater than that given by step 1 of the method 
Sieanalysis, the first value for deflection given in Tables El5, E18, 
and E19 does not represent the collapse deflection. 


For the partial loading conditions there is no direct relation 
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Figure 13. Load/Deflection Relation for the Arbitrary Plates of 
Tables El2 to El4 (Miid Steel Material Properties) 
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between the deflection at one value of the applied force and the next. 
Thus, once the uniform load is reached, the deflection should be within 
the same accuracy as the plates shown in Figures 6 and 7. This fact 
alone should help to substantiate the loading process. In all cases, 
the deflection/load relation, in general, follows a continuous curve. 
There are exceptions to this, but they are caused by the change over 
from the constant ice crushing pressure to the uniformly distributed 
moa If, for example, the deflection for the first partial load is 
greater than the deflections for subsequent loadings, the deflection 
will be decreasing as the total load increases. As can be seen in 
Figures 13 and 14, pene en leceen is never less than the previous value. 
Mimmeerore, if the accuracy of the calculated deflections for the uni- 
Pommily distributed loads can be accepted, then the calculated deflec- 
tions for the partial loadings should be of reasonable values for use 


iamedes failure analysis of the plates. 


~ 


Absorbed Energy Relation and Failure Criterion 


There is one basic purpose for calculating the energy absorbed 
by the plate during the loading operation. The calculated energy pro- 
vides a means of determining the possibility of plate failure on any 
increase in the load. There are several methods of determining the 
amount of energy absorbed by a plate. The method presented provides a 
Measure at the point of expected failure. 

It has been found in experimental tests of dynamically loaded 
clamped plates (see references E19 and E20) that a linear relation exists 


between the measure of applied force and the measure of the energy 
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absorbed by the plate. The relation between F and U in Figures 8 and 

9 is clearly linear for all plates. The results shown in Figure 10 for 
the arbitrary plates under the idealized ice-loads are not linear. The 
nonlinearity is due to the changing size of the bearing surface whereas 
the curves of Figures 8 and 9 are for a uniformly distributed load 
throughout the loading process. 

It is interesting to note that the only difference between 
the two curves in Figure 10 is the material properties of the plates. 
ieee versus U curve for each of the four plate thicknesses considered 
coincides with the other plates of the same material properties. This 
same Bepera 1 relation er to hold for the plates shown in Figures 
8 and 9, although B/A ratios and plate sizes differ in these figures. 
Thus, based on ae results of Figures 8 and 9, and the consistency of 
the results shown in Figure 10, the absorbed energy relation, U, should 


provide a good measure of the possible failure of a plate. 

The calculated values for’ 8 and ¢ from the data given in 
fepves ES to Ell are plotted in Figure ll. The values of U representing 
the maximum load that was applied to each of the experimental plates by 


(El) 


Loeser are also indicated. With the exception of plate number 
five, all of the maximum load points lie on or above the failure cri- 
terion line. Plate number five failed due to a weld defect, but the 
other plates did not fail in Loeser's experiments. Thus, there appears 
mempemet Least a rough correlation of the failure criterion for two 
different experimental results. 


Figure 12 shows the final application of the method of analysis. 


If the line of the criterion is accepted for the smaller values of § 
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Enem the thicker plates can be expected to fail shortly after their 
collapse deflection is reached. This is actually what would be 
expected for the very stiff plates. However, it seems that the line 
Of the failure criterion should include larger values of 9 for the 
smaller values of 6 . This is just a speculation, and the presently 
femened failure criterion anenee to provide a reasonable estimate of 


the failure load of a plate. 
CONCLUSIONS 


The method of analysis for large deflections of a clamped 
rectangular plate under idealized ice-loads provides a reasonable means 
of predicting the behavior of a plate, Furthermore, the solution 
method provides a simultaneous check on the possibility of plate fail- 
ure. The calculations are tedious, but not difficult, and the use of 
the computer program will greatly aid any extensive failure spiliershe. 

A failure criterion was Sought and obtained.  Furthemmore , 
the criterion is not limited to plates under ice-loads. It should be 


equally applicable to the analysis of a plate under any lateral load 


except in cases where shear effects are large. 
RECOMMENDATIONS 


The major recommendation is quite obvious, that is, to conduct 
ieee) LO provide more data for the failure criterion. Other areas that 
should be investigated are other possible solution methods for the 
load/deflection relations and the energy absorbed relation. These areas 
should be included in any failure tests since they are necessary for 


the establishment of the failure criterion. 
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Appendix A 


PLASTIC COLLAPSE LOAD 


The elastic moment at the center of a partially loaded simply 
supported rectangular plate is given by Timoshenko in reference Tl. The 
moment about the longitudinal axis at the mid-point is, 


M.. = 6P (Al) 


where P = qc and @ is given in Tables 20-22 of reference Tl. The 
values of @ depend on C/A and B/A, and those values given in Tables 
20-22 decrease in accuracy as C/A decreases. Values aoe G6 for the 
smaller values of C/A can be found from the equation for M given on 


tea of reference Tl. That equation is: 


ee CMC Ee) =) SL) Ca a-v)| (A2) 
a ek rc [5 


The parameters A _ and M. depend on B/A and the values for them are 
eee im Table 27 on p. 161 of reference Tl. The values of 6 paiem 
Wee from equation (A2) for C/A less than 0.20 are given in Table Al. 

In order to use equation (Al) for a plate with a clamped edge 
condition, the coefficient @ must be adjusted, since the values given 
in Tables 20-22 of reference Tl and Table Al are for a simply supported 
plate. The correction term @j1 may be found on p. 206 of reference Tl 
memevarious B/A ratios. Thus, the maximum moment at the center of a 
clamped rectangular plate under a partial load is, 


(Mo) max = (6+ 61) P= 6 P (A3) 


The factor 6° is plotted against C/A for three values of B/A 
ier teure Al. ‘Two steps were necessary in finding thesproper value 


of Qe’ . First, @ was plotted against C/A for the values of 6 given 





in Tables 20-22 of reference Tl and for those values of @ given in 
Table Al. Then for each B/A, a single curve was formed which combined 
the two solutions of @ . An example of the graphical solution 
for @ is shown in Figure A2. The second aeen mequrred adding the 
constant 6), to the new values of @ obtained graphically to complete 
the solution of G6" 

Once the elastic solution for the moment at the mid-point of 
a clamped rectangular plate was established, the plastic collapse load 
was found by simple plastic analysis. The plastic moment at a cross-~- 


S-eevon through the thickness of a rigid, perfectly plastic plate is 





My = To H? (A4) 
4 
Setting (4) max = M, the plastic collapse load becomes, 
Pennie | (A5) 
ao: 
eee Py, = gC? = PC? FOr a icereruchiae prea tne mequ baeonL 
c? me Op H? 
4 6° Po 
or Cs [oe (A6) 
2V6'P. 


By dividing through by the width of the plate, A, equation (A6) may 


be non-dimensionalized. 


a Oo “4e|(a)(2)f (A7) 
Be H A 


The desired relation between the load bearing surface and the ratio of 


material and ice properties for given plate geometry may be obtained from 
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equation (A7). 

The parameter O,/P, for values of B/A, C/A, and A/H has been 
tabulated in Tables A2, A3, and A4. From these values the curves of 
A/H versus O,/P, for various values of C/A have been drawn in Figures 
A3-A5. It was found that this was the best way to present the data in 
order to make use of them in the plate analysis. For a rectangular 
plate of given dimensions and yield stress being loaded by ice with a 
given crushing pressure, (C/A),, the size of the load bearing surface 
at the plastic collapse load may be found for A/H and O,/P, from the 


curves of Figures A3-A5. 
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Appendix B 
DERIVATION OF MEMBRANE FORMULA FOR PLATE DEFLECTION 


The solution for the deflection of a uniformly loaded mem- 
brane is given in several references. (See for example, references 
T1, T4, and El). The method presented here is the one given by 

ae 
Loeser. since the method has been developed elsewhere, only atue 
major points of the solution are given. 


The differential equation of a uniformly loaded rectangular 


membrane is: 


Ze 2: | 
Gee O2L q Bl 
Ox Oy 


where q is the lateral loading pressure and S is the membrane force. 
Wemmene plate dimensions given in Figure 2, the boundary conditions are: 


Z=OatxX=#0, X=A, y=rt 


Z = b, sin whe ae (B2) 


The b's are constant coefficients, and Y,s are dependent upon y only. 
After differentiating twice with respect to x and y and substituting 


for 97 and 022 , equation (Bl) becomes; 





Ox Oy 
Cog 
; Lee “ 
Demsun 0 fx Ne a ae ee eee (B3) 
Ra 2S _ 5 or 
n=1,3,5 - : 
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Expansion of the right hand side of equation (B3) into a 
Fourier series and using constant value for q between 0 © x § A 


results in the following expression: 


COD 


n-Ll 
q . 4a ; (1)? sin arx (BA) 
S tS n 1 
to ee 


Setting equations (B3) and (B4) equal to each other gives 





nol 
“7 Dre 2 
Yoes nett Veg giemng cr (B5) 
’ AZ ‘ S nib 


Solution of this ordinary differential equation for the conditions 


that Y, is symmetric about y = O and that Goa when y = + B/2 yields. 





n-l 
eee iG Ge cosh nitB _ cosh nifty (B6) 
Ores cosh n/tB 2A A 
Bet. 3b, 


2A 


~— we 


Then the membrane deflection may be obtained by substituting equation 


(B6) into equation (B2). 


oad 
n-l 
9 5 cosh nity 
Z = 4qgA_ (-1) A Sin Ditx (sb?) 
ioe. Ce l- -——--> = 
S wt cosh nfcB A 


2A 
Tf the constant terms are removed and the series is defined as k, the 


membrane deflection equation is: 


Z = lg? k , (B8) 


sn? 
and 


on 
a cosh nly 
k= 974355 (-1) hal ya Alok A (B9) 


6) A cosh nitB 
2A 


nN 


PD 





The values for k have been plotted in Figure Bl for three 


values of x/A. 


for the deflection. For example, the value 


W, of a uniformly loaded plate will require 
ing to the B/A ratio of the plate and x/A = 


feet hold for all values of q and S. 


The slope of the wembrane surface 


The use of these curves eases the burden of computation 


of the central deflection, 
the value of k correspond- 
Or DOF 


This same value of 


1S Obtained by diftteren— 


tiating equation (B8) with iespect to x. Thus, the slope at the mid- 
point of the long side of the plate is: 
CO 7 
: wo cosh nily 
OZ _ 4gé_ -1 Pesan x (B10) 
ox sft 2 n= 1,3,5 Ea cosh nivB A 
2A 
[Memory = O and x =O orxs 4, 
Co 
92 eee iat i nailed ae tee - (B11) 
ox Sire n? cosh nIvB 
2A 
Od 
Noting that ol “7 S . , equation (Bll) becomes 
Oo 
92,98 _ 4gd. 1 
= : iS 
Ox ZS sft: Ae L385 nZ2 cosh nitB 
2A 
or 
poe GA ka St x my =O (B12) 
Ox as 
where 


Ooo 
ki 4 e. 1 
d= Sa 


ye cosh nitB 


2A 


The coefficient k, is also iver in Figure Bl for values of B/A. 
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Ap endix. 


BIAXIAL STRAIN eLATTON 


Assume that a given element of the plate has only normal forces 


acting on each face. 





Sy 
ox m= Sx 
dz 
Sy* 
For the case of uniaxial loadfiip let § = § = 0. Then the 
Z 
strain components along each axis will bua 
CL = S Le eral = aay eae eal a -V 6.4 (Gil) 


where Cy eepresents the principal strauimeimehomeduccet Ton mine 


strain components in the y and z directiims are obtained by considering 


=S =0 dS =S = 0 respectively 1 
S : an : y PeSpecrive ya) Unuce 


x 
S52 —o © ee eee ee) 

and C63 = E 3; €,3 = =a, © 3: O,.3 = -VU €.3 (C3) 
Tf all the forces are acting, the total strain in the x direc- 


tion will be, 


lt 


€ 


x 


Cece ee Cee C | “VEL, al err 


if 


78 





For plates it is generally assumed that Cc. is small compared to the 
other strains. The total strain jn the x diveeuion for ¢ a O Shatlle jee 
Z 


x 


pam lar ly , tne total strain in (he y direction for oo = Owill Dee 


For the case where the in-plane strain is uniaxial in the x direction 


mee= 0 and, hence, €, =VE_ = a = 
€ y 9 al VE,. bLOmecaquaity on Cel) CL ec, 


ac, - Substituting into equation (C5) E = ( ee oe Ge, : CC 7) 
< 

This relation will be useful in obtaining the stress at a given point 

where the. strain in the x direclion is known, and the strain in the y 


direction is zero. The stress van be found from a stress-strain diagram 


for simple tension by using 


ie eae 
ier ae “- (C8) 


A more general relation is required when there is strain in 
pememetne x and y directions. IJ| ae “0, the strain in the x direction 


may be found by combining equations (C5) and (C6). 


oo ey “Ve ey “VE, 


I 


aa Va i oe 


S 


x 


(- v*) €1 -ve, (C9) 
When used in the form 
Gy = ex + VEy 
eA: ules 
the stress in the x direction may be found from a stress-strain diagram 


fem simple tension. 
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Appendix D 
PLATE DEFLECTION PROFILE AND AVERAGE TRANSVERSE STRAIN 


In order to use the membrane formula for the calculation of 
Weaee deflections, a relation between the deflection and the membrane 
force must be established. This is accomplished by assuming a plate 
Me@wection profile. Then, for a given deflection, the average strain 
feeoan the plate can be calculated. From the stress-strain curve for 
the average strain, a value for the in-plane or membrane stress can be 
found which is converted to the membrane force by multiplying by the 
peace thickness. 

For this analysis, it was desired to use an assumed deflec- 
tion profile that would represent the various loading conditions. 
Thus, the solution for a partially loaded simply supported rectangular 


(Tl 


plate as derived by Timoshenko was selected. The general expres- 
meeeeirom page 138 of reference Tl for the deflection at any point 


along the x axis (see Figure 2) is: 


Co m-1 
4, 2 : 
_ 4gA ()) sinm®@€ J) 4 i cosh(Q,,-2 7) 
ac Ts x DA cosh X_ 
sinh 2 @ 
cL @ psinh( KX 7 -2 a) + OY Ae Beg ace sin m/tx (D1) 
Zoos: Ce a 
es A 
m7TtB _ mitc 
where Xn 7 Ey and Oa OA 


Since only the shape defined by this expression is desired, the factors 


jaeiront of the summation sign were dropped. Then, in order to obtain 
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a non-dimensional profile, each value of deflection is divided by the 
deflection at x = A/2, the mid-point of the plate. The result is a 
series of shape coefficients which depend on B/A, C/A, and x/A. 

The calculations for w/W for various values of B/A, C/A, and 
x/A were made much easier with the aid of the digital computer. The 
short program used to make the calculations is included at the end of 
Mmieseappendix. Table D6 is a list of the input variables required 
for the program. This same program is incorporated into the overall 
ice/plate interaction program discussed in Appendix E. The output 
from the plate deflection profile and transverse strain computer pro- 
gram is given in Tables Dl to DS. A typical plot of w/W versus x/A 
is shown in Figure Dl. The variance in deflection profile for a con- 
Semerated load (C/A = 0.02) and a uniform load (C/A = 1.0) is illus- 
Mmeeeea in the figure. Other assumed deflection profiles are also 
feeetcd in order to erode a comparison with other membrane solutions. 

The only reason for assuming a deflected shape of the plate 
is to determine the average strain within the plane of the plate. This 
is accomplished by measuring the length of the deflection profile, 
subtracting the original length of plate, and dividing the difference 
Bemene Original length. If L equals the length of the transverse 
deflection profile, and A is the width of the plate, then the average 
transverse strain is: 


6, .-A. (in/in) (D2) 


In order to determine L from the w/W, x/A relations, the 


deflection profile was divided into a series of straight lines. This 
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O-EXPERIMENTAL DATA REF. €1 
B/N =2.0  q=1000 psi 
SEE TABLE Ell 


PARABOLIC PROFILE 







10 
CYLINDRICAL PROFILE 
W FROM 
W TABLE 09} C/A=1.0 
B/A=2.0 ) C/A=002 
0.5 
O 
05 0.4 0.3 Oe O.] 0.0 


xX 
A 


fmeeure Dl. Deflection Profile along the Transverse Axis of a Partially 
Loaded Simply Supported Plate 
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Iimicear approximation of the curve allowed rapid computation of L by 
the computer. In addition, Sag parameter W/A must be introduced so 
that the strain calculation is related to the actual Dice demlect Lon: 
Figure D2 represents a certain segment of the profile curve... By mul- 
tiplying each segment of A (w/W) by W/A, the are length AL will 


Eeemesent an actual segment of the deflection profile. 


=|z 


(W), 






DEFLECTION PROFILE 
CURVE 


A(w/W) 5 


Ww) 
W. -k I 


Figure D2. Element of Deflection Profile’ Curve 
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EAS , 


4 


A i = lA (e/a) | ‘ “- lA (w/W) (W/A) [ ()3)) 
where 

oom (BL (8, 
and 

A (x/A) = (i), - (i). 7 


Then, the actual length of the deflection profile will be: 
J 
L = ; i (D4) 
a ae J 
J 


where J equals the number of divisions of the x/A axis. With this 
value of L, the average transverse strain can be found from equation (D2). 
Por the case of a square plate ene sate race Gen ttt ahem, 
direction will equal the eae along the x axis. If B/A > 1.0, then 
the @vo strains will be different. The transverse strain will be the 
same as found from equation (D2). The longitudinal strain can be found 
using a similar relation, but the length of the cumyoumuce be adjusted 
to reflect the longer dimension of the plate. This can be accomplished 
simply by dividing the W/A term by B/A so that the parameter W/B is 
used in equation (D3) in place of W/A. Then, if Lt is the length of 


the arc along the y-axis, the average longitudinal strain is: 


Ey = L4 ae (D5) 
B 


Broa l ly sewer incipla straia ce seleemraopoimtwor Lherplaue 


can be determined from equation (C10) of Appendix C, which is repeated 
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fene £Or convenience, 
Ei 2 Ex "Vey : (D6 ) 
L oe | 

This is the value of strain that is presented in Tables Dl to D5 and 
Figure D3 for various values of B/A, C/A, and W/A. The strain is 
referred to as the average transverse strain in the tables because 
only deflections along the x-axis are computed in the analysis. 

Also, it should be noted that Poisson's Ratio is equal to 
0.5 for the calculations. This is the value that is commonly used 
Bemestrains in the plastic range. Both V = 0.3 and 0.5 were used in 
the calculations to test the validity of the membrane formula and 


V = 0.5 generally gave the best results. 


85 








fisrciet 


ets 
roe ats or CCE EEE eee w a tt 
. My Th ESE Seat Eppes Sea HH s(t STARE ee - Ht ate H ae PEE Sausesenedenane att ESSE Ag 
I r 


ieee 


: i HEH He erases | 
free ae EES Pe ee ! ce te 


eae a 





oe See aie SEER SEE CEEE CECE EE EEH et a SESH fete EEE ae sarees : 
eaem ah z Ht Ht Ph sWTHEE Pe tacky 
: natant ase a. 


et ae 
ae cE SLE PES Fey es SS 


Tae ae see - seuevesseeeune Perret a The sa esterase 
+ 


4 —t 
i 


ae Sine a HH Ede 
eee eer egrets 
Baas oa al ge 


a Bae eerera pat Heh satay efi He EH eee ee 


rie Pere Lh 
a SS es fee tba veut ea pee 
Pein Se setae ger eee ee ee eg ase 
ie S00 Pree ete 1 2 ieee 





2 ajo ge) ria etal ests mate ERueas aeaeeo Ee nee H sve a 
a L s HEH 


O eit ate See ee F CA ae fal ate | CES ce gace C1 Fete Pee 
| a ee SUT SR Ue ne 
BCE Pepe ee ae Be SHEE oo pote Cr 1p3 Hs HH H cE B E ene 


aS Ht 


font ae coo 
EES ESE EEE 
HHH Pt one LI ie Sint | ! a] | 
{+ 4 Sune aH eee ae HH AL a He -! ha ; 
EERE eae bp) Pp] fii + ci 
ae Hs ete + Ht PEER Ete Hh Ty | i HH a H a a poy HH ci aT 
oh ECE Ce rete oH EE oA EEEEEEEEE Eee oe CCEA Heo EEO Gscuecpeee 
1 Bw B || praetor Re ange iB. 4 


Ree EERE eee ott 
He oH : eae dena PEELE ecb 
SOR Paar a ar - aie AGE Pe cierrot 
ae) mt) 7 | Tf] cH to PE Baa eta e a ray ae t | 
HEE ane PCE HE: Ha eae a a eine oi HT HH Pty Soo Ni a ri He 
4 a Ht erie lalate ree Peete Be 


Eta Ae EE) eh ren gee Sm re 
o EE Eee caergos eae Ao Se Sait ee ELS 


86 





coo oH EEE eins : secnesesoanen H+ eee eee eee ee PH PEPE eee epic Poe BCH eee | eee 
ce HE es . — oo : eH sana aoe tse eae 
Ba PEE Eee FEE seas po + er <a Hart oo Hy} Lota soa 


imal sia ++ eke Ga 
Ba of a je aa aH a od fa aie 


op a et FE ro eeire EH SHH Sc See ate cary fj ees SS 
aoa e Soros He HEE Prete ft EEE a FF SERRE EEE EEE SSee cou ssevinetaces “serecesraseers 
Pooh apa te {a fp ae po ue SHIGE faa a a aaa rh ibis pera feet egabe ea fais 
Peet itt Eee eter al Eerste ES aa 7 FH Hee EEE ede | so] ea iui Sea aif peedsesacastoeee Tees rane : Saeteaiaie 
See Poe miciot taleleten giana er = ae BRezussenseceett POC tat fe Ti senenens 


ir fi a gv i | | 
Ga RSS RR See ici i eH HH HH eo ek 0 ge fm | oa | cw 
~ = oe we nami ale a al CeCe Pay le 
H a a a ne it ect totat= atcha Fa i a i Hi eV Tf ih fy i totatate Peet tt iat Fa) a gg eg Lg) 
inal a tian tamahea taal 


GSS 2 REREB EE eH EH fal] etetatal [age afc ea Ff kg eff A Tia a) |) ap 
& a Pt EES SR SS USSSA. SEER EEE EEE eee Bepee 
Pit tit Labatt atta tal aba batat tama =| = ~ a 
SIUERERA Fit maken famtat aad inl aa savant as] maT Os as Lan nT 
a ~— pooled ee ee eee Ane ne meine 





TRANSVERSE DEFLECTION PROFILE AND AVG, STRAIN 


B/A= 1.2900 DOT SS ONG Ae = Ot 


© 


wero OF DEF. COEFETC., At V/nemP ela Cmer LC fi Meer. 


Xa eo0 ao) S310 ees ao 
yA WiC MAX) 
OeorO 9, 00035 1.00@0@ 0,92356. @Q@572752 ©3504 4- Gis2> 
meee) 0.00141 1609000 69162 02.7278 6.5040 42549 
Meo! 0200316 11,9900 9.9171 0.7293 0,5949 05.2555 
Mme 0 0800559 1.0000 O491827 O.4.73212 060.5061 99,2564 
Meo 0.00862 1.9900 0.9194 1% 73325 0.5877 0.2574 
fee 05032325 wI000 0,928] 0.7497 6,51399 0.7642 
MmeO0 06115900 1.9099 90,9431 06,7828 96,4598 032228 
Ooo 0.20991 ##+«41499000 0.9505 9.8058 90,5807 00,3071 
MeO 0.783398 13,9909 O5.9537 0.48177 G.5999. 90,3177 
wee OO. 3107" 1490000 O08 9549 9.92713 6.6053 9.322739 
AVERAGE TRANSVERSE STRAIM FAR GIVEN MAX BDEFL. 

Me) 0 02500 N50) 600975) 1000 150) 
C/A 

mmeie oo) 0.00276 049011902 9, 0747C O484265 0,996443 2, 16232 
fee) 0,0027T4 9401103 0.9247? 0.04236S 0599677 9316847 
eee 9, 00277 C.01155 9.02474 0.04376 0.09699 9.16869 
mee) 0, 00277 0.91107 3.907480 0,04224 0,9097902 0,14692° 
cree, 0,00278 0.61110 0697486 0504395 0.997350 9514935 
mmm 4. 00223 9,01178 0497578 6.94467 9.099°33 0.171846 
eee) 0.09225 O04 C1175 960263? 0604647 9.10244 9.17208 
fee 0.00395 2,01217 7.02723 6.04804 0,19591 0419335 
mm 00312 0.91245 0692725 O.04°11 Of19810 9.13425 
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Table D6 
LIST OF INPUT VARIABLES FOR COMPUTER PROGRAM FOR PLATE DEFLECTLON 


PROFTLE AND TRANSVERSE STRAIN, 


Computer Definition, on oymoo (in Text 
Input Variables” 
BAL First value of B/A to be used in Eqn. (D1) 


G2 Increment of B/A to be added to previous 
value for next set of calculations 


63 Last value of B/A desired to be used in 
calculations 

PR Poisson's Ratio 

NWA . Number of W/A values to be used in calculations 

NXA Number of x/A values to be displayed in output 

WA W/A 


af. .f. 
ay 


Output Variables” 


W (MAX) Solution of Gamat i1one (Dl) stor x). -a0 5. 
the mid-point of the plate 

x/A For plate qeneiee Ch Le ure ke 

C/A For load definition of Figure 1 

W/A hor defleculonmdertinicien or hi etnecms 


of, 
a 


See statements DEFLO0O04 and DEFLOOO6 in the program listing. 


See Tables DI - D5. 
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Appendix E 
COMPUTER PROGRAM FOR ICE/PLATE INTERACTION 


The FORTRAN IV computer program to calculate the mid-point 
deflection of a plate and the amount of energy absorbed at the expected 
point of failure is listed at the end of this appendix. The logic of 
the program follows the step-by-step procedure given on pages 37 to 39. 
The values obtained for (C/A), from step 1 is used as input to the 
program. The applied force F, and its corresponding W are determined 
ireeae Lirst part of the program, Then the absorbed energy U is cal- 
Gamated from the esas relation at the mid-point of the long side of 
Mme plate. Next, the applied force is increased by a AF and the pro- 
cess is repeated to find the new value of W and U. The output format 
Geeeme program is illustrated in Tables El to E19, The description 
@emene input and output variables and their relation ee the parameters 
Seiied on the page of Nomenclature are given in Table E20. 

Similar to a handvcalcul@tionsof W and U, a cemstitutive 
relation between the stress and the strain for the plate material is 
meguered for the computer program. A generalized stress-strain curve 


was used in the program to obtain a stress for a given total strain. 


The general form of the equation is: 
O=vDC A + €))) = De), (El) 


where € 5 ig the plastic Strain and) © -icstrwectal strain of the 
Memeerial. With the proper choice of coefficients, this equation pro- 


vides a reasonable fit to an actual tensile test stress-strain curve. 
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It was found that by setting A equal to the yield strain of a given 
material, and using the boundary condition O = on when € = € the 
exponent n and strength Coefficient D can Be determined, For materials 
such as steels which show near constant stress for the smaller values 
of € , the value used for X should be the strain at the point where 

OQ ceases to be constant. For example, the stress-strain curve for 
mild steel exhibits a near constant stress at the yield point up to 
approximately € = 0.02 in/in. Thus A = 0.02 should be used in 
moamarion (£1). 

Since the generalized expression is always increasing with 
increasing Sa tie Wehoiecemar ce . is important to ensure a proper 
curve fit. It is usually difficult to determine from an actual stress- 
strain curve, the actual value of strain that corresponds to O,,. 
Metemeest selection of €,, is at the point where the 0 -€ curve first 
starts to level off. The slope of the curve may not actually ere. 
me macadl just prior to the failure strain. But the ultimate stress 
is usually approached to within a very small error at a strain con- 
siderably smaller than the actual Cu Therefore, in order to main- 
tain a good fit of the curve between © = A and € = or a value 
of €,, should be used that may be less than the actual ultimate strain. 


From the chosen values for A and € y> the values for D 


and n are determined from: 


Lit: ¢./ 1a (E2) 
= eee 
rn Gey eee) 
and 
Gi | 
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The use of the generalized stress-strain curve allowed the 
incorporation Be ithe average transverse strain program discussed in 
Appendix D into the overall program as a subroutine. The name of this 
subroutine is STRAIN and its listing follows that of the main program 
at the end of this appendix. For an input of the variables B/A, C/A, 
and W/A to STRAIN, a value of €1 as given by equation (C10) is cal- 
culated and returned to the main program. Then the stress for € 


is found from the generalized stress-strain function. 
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ICE/PLATE INTERACTION 
ICE CRUSHING PRESSURE= 45,2 PSI 
PLATE DIMENSIONS (INCHES) — 
TeMeKtiEss— ce 


WIDTH= 24,00 
LENGTH= 84,00 


MATER TAL ROP Eig bes 

YIELD Sie SS= 5600050 PSa 

UL. STRESS= ) (22a ew 

YIELD STRAIN=9,0200 IN/IN 

LT,  STRATSeOs 1400 IN71* 

POISSON S RATI]U=0,500 
AT e-PRO LONG STOE 
i W S C/A STRAIN STRESS U 

Hee S } { IN) (LR/IN) (CIN/SEN) (PST) Cie 
Meeeee?’s: 11,9068 T0000 1500050058 2305 >> -Giy> 1309.28 
eee i,) 3,39133 #O000,0. 1,000 0,0460 62054. Ze ae 
memoe?,) 5,2839 T5717,.6 1,090  Qe@064527 64903, 2 47 FO 55 
Pee ?,) 68,652 RO1S.4 le O00 Os Oe. 9/41.,0 


Mmeeeeo2’,i f.9763 836652 1a 3OOCNO yy 242] ees ZO eee 


Tablé El. Plate IB of Reference E3 
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(LBS) 


melo 10.3 
010.3 
ee25810,3 
Pero 10.3 


11€S81C.° 


ICE/PLATE INTERACTION 


16 CRUSHING 


PCATE DIME NSTONS Sy cies. 
TT CRNES S = erie e 


WY OT H= 


LENGTH= 


MATERIAL PROPERTIES 


YIELO Sisesi= 
SEs Si 
YIELD SPRL oe oie S 


ULT. 


ULT. 
Pols oie 
W 5 Cua 
(IN) CLB/IN) 
Pe39 8) 1551567 12000 
Bis £318 CT42.2 See 
4.5€42 BOIG6 Dale ae 
Onc) 9437 6S ieee 
Wel39G  10049.53 Wienoan 
Table E2. 


S257 0% Ohi | 


IN/IN 
STRAIN=0.1600 IN/IN 
RAT I9=9 2509 


PRESSURE= 48.9 PSI 
Sao 
S400 
41800.0 PSI 


AT 41D PT. OF tN Case E 


See 


(IN/IN) 


One 33 


0.0660 


Ors 933 


O. 1194 


0.1449 


vie 
ist) 


47631.6 


54727.4 


DSi Oe 6 


64197.4 


Plate 2A of Reference E3 
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ICE/PLATE INTERACTION 
ICE CRUSHING PRESSURE= 241.0 PSI 


PLATE DIMENSIONS (INCHES) 
THICKNESS =0,244 
ITOTH= 13,50 
CENGTH= 2c 


MATERTAL PROPERTIES 
YIELD STRESS= 54000,0 PSI 
LT, STRESS= 70700.0 PSI 
YIELD STRAIN=0,0299 IN/IN 
ULT, STRAIN=0,1400 IN/IN 
PIUISSON S RATIO=0,500 


Ai orp 2 eee ee Coste 
= Vf S C/A Gear sues s U 
meas) (IM) (LB/IN) (CIN/IN) (PST) (ios 


Meee oo? 0.232755 13664.0 1.000 90306 53923. 180 5a00 


IND 


Meo. 0.6505 1346460 18000 0.9612 64028,4 Bao oer 
emer ss > «€©60.9714 1372081 %12.909 0.0916 47194,28 Ola ayo 
fee so) 86122257 1449661 1,000 0.1186 693C7,7 821352 


mee od ln 4683 6 615125,1 12000 0.3450 7O997,;8 10294,8 


Table E3. Plate P2 of Reference E2 
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[CE/PEARE INT Iah 


ICE CRUSHING PRESSURE= 19020 PSI 
PLATE OIMENSIGNS ( TMC hes} 
THICKNESS=021)49 
WIOTH= L320 
LENGTH= 21,00 


MATERTAL 2 RGPEr ii. 
YIELD STRESS= 56000600 PSI 
ULT, STRESS= 70ers 
YIELD STRAIN=050799 17) 
ULT. STRAIN=O,1400 IN/IN 
PJTSSOM S RATT O=Osee 


AY “ID=P7T OF LOMIG SIDE 


E W S C/A SN SSieess 
CER S ) (IN) (LB/IN) (IN/IN) (est) 
meee, 0.526? 566450 Us000 95929) S57273-4 

memes, 0 15,9332 678359 laos O CAs 62 Vous 7 
/meeo>,) 11,4329 7335,0 7 13000 903 074325 36402 
emeoo,]) 142809) (74654 Is@QO OeiGees Geeo0.1 
Meo.) 2.1707 8069.4 1,000 “Qeelnos “6262554 
Tepe ES. Plate P3 of Reference E2 
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Uj 
(IN-LB) 


133 2a! 
C2 ieee 
4469, 3 
522 Ges 


Tee x 





E 
(LARS) 


34070,0 
68070,0 
mee 20.9 
poee2zC.O 


Pewee O.0 


LCE Pie Tee 7 ON 


ICF CRUSHING PRESSURE= 32020 PST 
PLATE DIMENSTONS (CINCHES) 
TIT Cai eS S- Gas 
WTS 13500 
LENGTH= 21.00 
MATERIAL PROPERTIES 
VIEUD SIPeSS= 47060205) 4) 
ULT, STRESS= “4GHe@GwO =? so 
YIELD STRAIN=0-0200 [17 3) 
ULT. STRAIN=0,1600 IN/IN 
POISSON S$ BAT IG=Ge500 
AY MYID-PT OF LONG SIDE 
W “ $s CaN Sue nIN STRESS U 
(IN) (LB/IN) (IN/IN) CPS} { IN-L®) 
OO 7? 3360,0 O00 Oe 0267 wees oer. f B9R2,4 
ire 20 265356 TeO0 me Oe! 37 toa 1925.3 3 
ieee. 405458 3000 iO. wa oreoe 2 ZOLISZ 
Peat) 7 4360,6 WIOOC MGs C33 Gua 257 392653 
Pea) 4671,8 We8ooo 96,40 ease eu 4950,?2 
Table E5. Plate P4 of Reference E2 
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TCE/PLATE TNIERACT ih 
1C™ CRUSHING. PRESSURE — Bie0 PST 


PILATE DUMENSTONS (INGE se 
THICKNESS=0,068 
WIDTH= yer oo 
LENGTH= 21,009 


MATFRIAL PROPERTIES 
YIELD SE@ESS= 22000).0 oa 
ULT, STRESS= 418060,;0 PST 
YIELD S§BAIN=Omb200 IN/ITN 
ULT, STRAIN=0,1600 IN/IN 
POISSON S RATIOQ=0,500 


AVS Oe eo Veh “Saat 


fe W S Ci AN Sait eS Gs Ss U 
GieRS) (IN) (Liaise .. (IN/IN) eat (IN-LA) 
Meet >> 00,7379 2171630 ASG OO SOS wee Gas 616, 6 


ees > =, 3536 2370,4 1,000 030370 34627,5 1280,9 


M0g4,5 1.8799 259755 V0 0897 oe) 193 ee3 


“J 
N 
poms 
Wh 


pees > = =—2 3545 27ga@4s,2 \rO0OO OTe 2 4 > 4 2585,.8 


fees ,5 «62,8169 2846,9 159900 00,0844 38502,4 224956 


Table E6. Plate P6 of Reference F2 
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YOE/PRATE Vier Ae] 1 


TCE CRUSHING PRESSURE=  224,2 es] 
PLATE DIMENSTONS CINCHES) 
THICKNESS=0,119 
wrOTH= 12550 


MATERIAL PROPERTIES 
YIFLO STRFESS= 42000.0 PSI 
ULT. STRESS= 61000.0 PSI 
YIELD STRAIN=020200 IN/IN 
ULT>s STRAIN=0,1600 IN/IN 
POISSON §S RATIO=0.,500 


eee ee Sree 


F W S Ca STR ©Sees U 
Cie S } (1M) CER? in CIN/IN) Ces) (iN=28) 
mee, / 0,6802 4998,0 1,000 040312 45492, 8 1420,0 

104420,7 11,2409 5478.4 1.000 060595 51074uani S03 is 9 

meee CO, f 1.6747 6088.3 1,00G 00,0846 S44] 0e? 4604;6 

Meee O,f 230776 655Gb ITE OOCOTN Gn ROe° 5627557 6198,° 

epee, f 2,444? 6952.4 bar Oana 5a es TBs 2 
Tapmepn/, Pilate Ps of Reference E2 
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LGG7 Pie Seer Gr ay 
ICE CRUSHING PRESSU?E—- = 100. oe 


Pigott “DUNS FOMS ste Grae) 
THICKNESS=0, 188 
WIDTH= 13,00 
LENGiIH= 2a300 


WATEM LAL! PR@OER WES 
Ve ee oo Oe 
ULT. STRFSS= 5900020 PS! 
YIELD STRAIN=0,0200 IN/IN 
UIT. STRAIN=0,1600 IN/IN 
POVSSON SS RAE O5 es 


AT i= Pied LONG Soe 


F Vy) S (fn STRAIN STRESS UJ 
(LBS) (IN) (LR/IN) (IN/TN) (PST) (IN-LB) 
464800,0 9.2842 7238,0 1,900 0,0205 38691,9 192 of 


MaaoOeC 068524 7238.0 1.000 0.0614 48472,2 2977.2 
223800.9 1.2766 8953.7 1.000 0.0968 532146,3 5151.0 
327390.0 1.6229 §8869,0 1.000 031292 56462;5 729245 


420800,.0 14°440 901989 TeOO00 OF GOS” 2904 leas 947941 


Table E8. Plate no. 2 of Reference El 


- 
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F 
vies } 


25625,9 
58825,0 
92625,0 
126425,C 
1602725,0 
194925,0 


eee 5.0 


IND 


Bono 25,0 


7 5.0 


LCE rer 


PLATE 


eee 


i] 
(IN) 


Ge 2265 


ORS SY oy 5a 


OS 252 


is ea 


Mie 0 Ree 


ero OO 1 


1,682? 


ia 


N 
> 

Oo 
aN 
—) 
oO 


ICE/Pi are. Pepe eer ON 


LOG > Oe! 


UST ING PRESsu 


DIMENSTONS (INCHES) 
THICKNESS = Gee! Gat 


WIOTH= 13200) 
LFNGTH= 19,25 
Ae lees 


YIELD STRESS= 3850009 PST 
WET. STRESS] Sc. et 


VIELD STPAIN-Ge070) Sila) 
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ULTs STRAIN=0,1690 IN/IN 
POTSSON § RATIN=0, 509 
Ka aes OF {ONG SUDE 
eS Coa SPRAIN 25 1 Ress lJ 
(LB/ IN) (JN/IN) Ces) Cee 
7232220  TaOO Os) i ee a5 Oe So 
7238-0 WeOOr On ec. 18/4,6 
7238,0 TGR O> 30. Coos Fao 1: Oe 
7 Re 1,000 0.0864 51983,4 44nn-—% 
937126 Ts000 BOs oo) face ee eee 
PALS h  i@O0 MOane2 16 26044. 7 691, 3 
9942.3. G2 900 eee 57660, 8 B249,5 
9635.4 “Je Gt Ge oe Oo. Soe 8 
997421 pene, 0st (92 655 Sir ones 
Table E9. Plate no. 3 of Refeiinee El 





(C27 Pile. Gee eC 1 ie 


[Ce CRUSING Pibsoun 


ti 


= Boee 2 P54 


PLATE ERNST ONS (hea s.) 
Ti TORN SS ye 
WY OTH= 1 Seger’ 
LENG H= Titeuhrc 


PATEK a ree ee 


YIELD STRESS= 23502350 @25 1 
Utes STRESS= GSO eee! 
YIELD SRAM Oo Oe. 
To SI Ole ey 


Pel Sot) Se a lee eae) 


AT MID-PT OF LANG SIDE 
E WW S C/A Ce eee S U 


(LBS) CIN)  (LB/IN) CIN/IN) (PST) CTN-LR) 
Meeo.0 0.1717 723819 1.002 0.0142 3587¢.3 508 2 
Mo.) 3.5151 7238.C 16009 0.0425 4464459 1910,5 
Mec.c 4600.8482 460s 7325.9 1wO0O OFOTO4 49@4R.0 5ie8.¢ 
Meso.) «21.074 8031.4 16999 9.9935 52955.1 4949.2 
Memiec. 8 «61.2245 3490.0 Pep00 9.9962 55247.1 6417.¢ 
MeeOc.0 (1. 4866)6«6 $155.1 «16986 0.1382 $7241.9 7903.8 
Mes.t «1.4752 9545.2 Jy290 0.1566 5896.2 $408.23 


fae. C deen) 0 Lee ee aey Oe ee. eo) 7 e518 ra @ 


Meee. «862, °2596 1f40S,6C 1.00C “O,2015 Gilbs27.€ 12407.5 


Table E10. Plate no. 5 of Reference El 
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. 
Ces ) 


33800,0 
1013000 
168700,0 
236300,0 


Bo 8 OQ 


a 


371200, 


oe) 


ICE7PUNGE Teac Tt iloN 


TCE CRUSHING PRESSUME] See. o-Ps! 


PLA@E DIMEMSIONS (JMEHEFS) 


THIEKNESS=@e 250 
WIOTH= 1 39100 
LENGTH= 26,00 


MATERIAL PROPFRITTIES 


O,1167 


Of 


WJ 
(TN) 


Oso er 


0.8157 


lee O28 9 


1, 


28TS 


VYIEDDE STR ss= coeur eee. 
UIT. STRFSS= B8360Gm0 PST 
YIELD STRAIN=90,90200 IN/IN 
ULT, STRAIN=0,1200 IN/IN 
POTSSON S RATIO=0,500 


AT MID-PT GF LGNG SIDE 


Ss C/A STicath) See ss U 
(LB/IN) (IN/TN) Cen) Cia) 
1650:0,0 ol COCR. 0! bo tere Z. 69954 


16500,0 Vonooemos4n rOsco.7) 241s 
16500,0 1.000 0.0569 75439,6 4290.5 
14500,0 Wage Osemes Gemeien4 Goon 
Weolme) 1. 000ml (Cl sumntegano eames 


17513,8 Is@00 “Oo leee2 8232107, 2 Wowmer ) 


mabue tel. Platesno. 4 of Reference El 
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Ee 
(LPS) 


423232600 


47325.9 


ye 3355° 


5.9 


4472325, ° 


35.9 


642335.9 


ICE CRUBBING PRE GSUar= 


Platte 


MATERT 


W 
(TN) 


O02 8001 


La°7erl 3 


1, 3359 


1. S556 


se 5424 


le 7079 


Pele 


IGE/P CATS EN GemecTi On 


DURLENS OMS. (1 Cais) 
THICKNESS S=05750 
WIOTH= 24,00 
LENGTH= 48,00 


Al PROPERTIES 

YTELD SYRESS= 2250020 Ps2 
ULT>s STRESS=» SS@ge. oO. P sy 
YIELD STPAIN=0.,0200 IN/IN 
ULTs. STRATHNRE Ow COO LI N/a 
POISSON SRAM =0, 520 


ATS ID=Ct Or aie Side 


S C/A STRAIN STRESS oy 
(LA/IN) (IN/IN) (PST) (IN-LB) 
28875,0 0,350 0.0757 50595,7 382944 
TRG SulO 05642 MOR ee Uemeols5 weeps 
278875,0 92837 0.1811 69518,4 1095832 
P8A75,0 ONSG5 Ou2libe Semee.s 124905—e 
2287520 16GD0 MOM 46 Poco) 15 acme 
7BR15,0 BeaGOO “eee Veseene 0 iaeagee 
29316,0 WACOM Os 2-ummac 7 5. n7 no oeee 
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Dapwewim2?. Arbitrary Plate no. l 


PeSe BATE INTee ACTION 
ICE CRUSHES SsUye— 7600.0 PS! 
PLATE DISENSIONS (iver s ) 
THICKNESS=1,000 


WIDTH= 24,00 
LENGTH= 48,90 


MATERIAL PROPERTIFS 
YIFLD STRESS= 38500,0 ST 
ULT. STRESS= 59000,0 PS] 
YIFLD STRAIN=0.C200 IM/IN 
ULT, STRAIN=0,1690 IN/IN 
PIISSON S RATIO=0,590 
AT MID-PT OF LONG SIDE 
F W wei C/A STRAIN STRESS OU 
CLRS) (IN) (LB/TN) (IN/IN) (PST) (IN-LB) 
104543.9 00,8259 38590,0 9,550 0.1189 55515,0 66092,8 
204543,9 0.9756 38500,0 0,769 0.1664 59474,3 9894,4 
3N4543,9 1,021] 38500,9 03939 9.2036 61954,5 12576,6 
494542,9 161063 39500,0 1,009 0,2340 63786,7 1492378 
504543,9 1.2355 38500,9 12900 0:2613 65249,.5 1704828 


eee? «686i 35240 «6©638500,9 1,900 O52260 66471,7 19016 


Table Hi3. Arbitrary Plate no. 2 
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fs 
(LRS) 


279935,7 


479925.,7 


479935,7 


Bee), / 


eee 5, | 


ree CE 


PIR 


MATER] 


W 
(IN) 


a) 


Onieni2 c 
0.8577 
90,9640 
| Geter 


11,1474 


ICE/PLATE INTERACTION 


USHING PRESSURE= 600.0 PSI 
DIMENSIONS (INCHES) 
THICKNES S= 15 250 

WIOTH= 24,00 

LENGTH= 48,09 


A PR@P ER TAS 

WIERD STRESS = SB5005.0 251 
Utd, STRESS=  .59070,0 PSI 
YIELN STRAIN=0;3;0200 IN/IN 


ULT. STRAIN=0,1600 IN/IN 
POTSSON S RATI9=0,500 
VISCO SPI AOr NG: soe 
. Cik ST Rea STRESS | U 
(LB/ TN) CIN/IN) (PST) (ri 8) 
Oxt2oe0 05000 041946 Gi421,0 17954520 
SGN2550 \e@OO0 O.2257 63377,1 426080 
4£6@125,0 1,900 0.22548 64915,5 16542,7 
42425,0 14.000 95.2801 66158.8 18541,4 
Geveees OO OC) O54 2035 “Giciie>, 20406,¢ 
Table E14. Arbitrary Plate no. 3 
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c 
(LBS) 


eee 00.0 
pee 00.0 
#91 200.0 


09,0 


Leer er 


PEATE 


hair ea 


W 
(IN) 


OFasT 7 
12 eer4 
1, CS4a7 


1 » i: 545 


ICE/PLATE INTERACTION 
USHING PRESSURE= 600.0 PSI 
DIMENSIONS CINCHES) 
THICKNFSS=1,500 


WIDTH= 24400 
LENGTH= 486,09 


PE CROP er iene 

VIELE SSG eas s- 320.05 s! 

its “STE Ss- ~Seeces cars! 

YIELD STRAIN=9.,90200 IN/IN 

ULT. STRAIN=02.1600 IN/IN 

POL SSON GS Rael). 500 

Wie Vari eo LONG Shoe 
S CTA STEIN STRESS UJ 

(LB/ IN) (CIN/IN) (PST) ( TNSwe) 
Sexo, Ot SGO0 O22162 6276450 24577356 
Dimes leer 0, 3272  “G833442. 2235568 5 
Ses 0 12000 05.3473 4917221 2mOZ i.) 


Diente OO Vee 8064662 699272,4 2561129 


Table E15. Arbitrary Plate no. 4 


ar 





[CEPR ES Tera iON 


ICE CRUSHING PRESSES —- GOClOers! 
PLATE DIRMENS BONS  (iNGres? 
THC eS S= 0077 30 
WIDTH= 24,00 
LENGTH= 48.00 


MATER TASC eer Riess 
YIELD SiSESs= 660000 0S) 
ULT, STRESS= S8300050 Ps! 
YIELD STRAITN=0,0200 IN/IEN 
ULT, STRAIN=0,1290 IN/IN 
PUTSSON S RATTO=0., 500 


Kee NG Se 


F Ww s Csi STRAINS STRESS U 
Cla } (TN) (LR/TN) (IN/TN) (ost Cale 
meee sooS 0.6423 49500,0 0,550 0.0694 77382; 3 936655 
eee? )«€©6C0L 783306495900,0 0-858 0.1083 81914,0 SiO o 
ees?) =€©600.8605 «849500,0 1,000 931365 84377,.4% 31515,8 
B989435°9 1,0074 49500,0 1.000 0,1598 860°66,6 13757354 
Pesos 9 «6141355 49500,0 1,000 00,1801 &7425,9 15746,1 
fees?) «61,2506 49500,0 14090 0,1984 88510,9 17556,9 
Taplesio. Arbitrary Plate no. 5 
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POE, VATE SINT GRACT ION 
ICE CRUSHING S Rees URe= ~GOO.0 Ps I 


PLATE CIMENSIONS (INCHES) 
THIGCKNES S= ies Oo 


WIDTH= 24.00 
LENGTH= 48-09 
PA Ro Nees B® seit ES 


YEO GheSs— G6000 20 seal 
WE Ties poepes S= Sea 0.9 O gi 5] 
YIELO STRAIN=0.0200 IN/IN 
ULT. STRAIN=0.1290 IN/IN 
POS SONS ai ne) tere) 


Ae D=P 7 Or WON Ges Oe 
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i W S Cy Sia Neo BRE SS U 
neo 5 ) (IN) evel) (INZIN) est) GaN UG) 
mel63.9 O65765 66000.9 08800 Owl009 8118!1.3 GTS 205 
Peeecse 9 O46182 G660C0..0 412000 O24 1307 $391%.2 IJCSIC.2 
Meese S O8F.1325 «6©66000.0 1,000 O41549 .65755.7 13284.4 
671182.9 0248313 66000.5 1.000 0.1758 87154.2 15321.2 
Se 8 he69 «6069195 66900.,0 1.000 O401944 88285.8 17167.C 
peeetciwo «=6©60.9999 66900.9 1.090 9.2115 89238.2 1838370.6 
Tapa 7. Arbitrary Plate no. 6 


lCE/PEAT I NaER ACTION 
ICE CRUSHENG PRESsUer- 1 600,00 PSI 
PLATE BlVEN Sa ees) 
THIGRNES Sade 250 


WIOQTH= 24,00 
LENGTH= 48.00 


MATERIAL PROPERTIES 

YIFLD STRESS= 64090,9 PSI 

ULT, STRESS= 8300020 PSI 

YIELM STRAIN=0,0200 IN/IN 

ULT, STRAIN=0.1200 IN/IN 

POTSSON S RATIO=0,500 
AT MIO-PT OF LONG SIDE 
F WW as C/A  STPAIN STRESS U 

(LES) (IN) (LR/IN) CIN/IM) = (PST) (IN-L8) 
200.0 0.4772 82500,0 1.000 0.1261 €3531.9 10537,1 
891200.0 0.7663 82500.9 1,900 0.2075 88749,0 17977.7 
POW1200.C 9,8479 8$2500,0 1.900 0.2241 99905,7 20152.2 
1291200.0 9,9724 82500,0 1,000 0.2438 90878,6 2215824 


Table E18. Arbitrary Plate no. 7 
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F 
Vee S ) 


ere), 0 
e277 00.0 
ipo 2 00,0 
meg 1200.0 


fo 7 00.9 


EGE Cs 


PLATE 


MATER] 


W 
(Tis) 


Oe er 
O32 56 
OQ. (066 
0, 7686 


Ora aCe 


LCE7PEATS ENTER ACTION 
USHING "PRE SoURE= "GO0s07P5 1 
OIPMENS Peas  CriCuES;? 
PGPEKNess—15 500 


LENGTR= 48,00 


AL PROPERTIES 

VUELOS So] 66 000, od 

Ulloa SIRES =) 830 CORO a Ss 

YIELD STRAIN=9050200 IN/IN 

ULT, STRAIN=9,1200 IN/IN 

POTSSGN S RATI9=0., 500 

Af MID-PY OF LONG SIDE 
S fe STRAIN SS Re Ss U 

7 Te (IN/IN) (PS) { IN-LR) 
SI00050— 1Ls000, “Onl261 835232159 10534. 
COO OO ses OOCNN Gm2025 8874950 1797757 
$9000,90 OOo Os27 1 892905, 7 ZO AES 2 oer 
SIVOC Oo ts UCC 04 2438 5 906738. 6 2215Res 


SOOO. Ve OO0 mms ci 2) PIOZTGe st 7508254, 


Table E19. Arbitrary Plate no. 8 


ey 





Table E20 
LIST OF INPUT VARIABLES FOR THE COMPUTER PROGRAM 


FOR PLATE DEFLECTION/LOAD RELATION 


Computer Definition or Symbol in Text 


+ 


Input epee eee 


ELD 0, 

ULT oe 

EY i A in equation (El) 

BF a 

PR Poisson's Ratio 

DF vet 

FLIM Last value of F to be used in the 
calculations 

H H 

A A : 

B B 

CA (ci. 

PC Ee 

Q k 

QD kg 


Output Variables (See Tables El to E19) 
All ef the variables displayed in the output are self explan- 
Peemyeon OL Wehne same notation as used in the text. The STRAIN at the 


mid-point of the long side is equal to ©, in equation (19) of the 
texmeee ine STRESS at the same point is equal to O.. 


See statements FATLOO03 and FAILO005 in the program listing. 


3} 





9c0U0 TIVGA 
Se d0TIVA 
yYVEVITIVG 
EcOQOTIVGA 
Geol tive 
Teool Vs 
JCOOTI VS 
6200 TVA 
SZOOCTIVS 
£200 TIV4A 
92001] Va 
Sé00 Viva 
COU yes 
CcOUTLVG 
ccOo0rlVvd 
LevotvIVa 
OCOe mea 
6LOOTVG 
BTOUTIVA 
L£TOOTIVG 
9100711 V4 
STOO TIVG 
STOOTIVG 
SloCFl va 
CLOGTIVaA 
tloollve 
VtU0G TI Va 
SOU Gta yS 
8QOU0 TIVA 
£GO001TIV4A 
SOC 2 
Scoellva 
YOUO 1IVG 
EQ0U TIVA 
cOoO TI VSG 
10GGiive 


(Sxeld) /(Wd%*O0% ( CK«{ O° 2/9) ) xO° ST) HWZ 


"Xo “3 (1S2) 
“XU1/ * sie ne Goo dts 


OLT {Ott ao 


O° T=43 10H 
O° O=XIC1 I 
H«xO V3 TA=S 
3d=1d 

toe leOd=Z 
J= TD 

Gil ages 
G=19 
(9°T-V3) GI 
VxV9=9 
V/G=Va 


(dM dx*Ad )/G T3LA=090S 
(453 /ARZ) 90 1V/ (1 107 GTA ieee 


(NI/NI)e*XS *s (NI/G1) 0 *XE4s END} 
NIVHLS 0 °XE a7 oe ee 


tcwy~OT7 T° c=eld 
(a(dI-NI}. € 


r*° XS ORCS Gane ec 


“ams *X8S* ade’ XeT/*edd1S YNOT BO Ld-GIW LVs*XS> Jey ivr 
(o2e€°) 3 ll ome 


"KOC “ais NE gia > So ee ce 


Pa=NiVdlS GialAe*X9?/ “HIS d el 6d asa es 


*T°8d*e=SSSYLS GISTAs *X92/* sS3TLeddOUd IWIYSLVW. 
Ud>a4 74 2G eee. 
(2 °2L5* e=HLONA 1s °X92/ *C° 8a * a= 0 ee 


*X92/'C°S3 Se =SSINNDE AL eo * KOC/* » (SAHIN) 


(C°os* ,=OI1LVeY S NOSS] Oda 
"EW * X9C/ 2 ae NN 
"L1Nse*X9C/7 el Sd tee 


yea a 


*XOc//ILVWUOF 


{ 


GTE’O} I Lig 


SNOISN3W10 Jd1LVi dae 


*XOC//*.ISd « ST ° Ld %2 =SUNSSSdd ONIASNUS 49 
*XOC//*F ANOTLIWGSINI SLVid/3312°XO0C/// “thls v7 ae 
(SOCf9)SLTUM 
(GO Omer Wage 
@® * O* Dd* V9 a* ve 


G*y *Ht od 


GLE°G0S ‘00S 


NIV d* 30° Ud*® 4a° Ads Ga. 


(Ove *s)avar 
(QIZ1A) GI 
(37 §S )OVAd 
(UMd xx LS) *xODS=CLS)AVWOIS 


NOILV isd OWO l/ NOT ageRe ieee id 


OQEe 


STtt 
SOT 


OTT 


Ge 


GO 


7c 
GLE 
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cLOOIIVG 
PLO 1 Va 
OLOOTIVG 
6900 11V4G 
BIODOTIVG 
£980TI VA 
9900 TI VA 
G900 TIVE 
YVIODOTIVA 
£900 11V4 
CYO00TIVGA 
L9U0TI V4 
O90ITIVSa 
6S00T11V4 
8SOO0TIVA 
£GsOu0TIV4 
9S001IV4 
GGUOTIVG 
4SOOTIVG 
esvoliva 

svO1iVd 
TSU0TIVGA 
OSOUTIVA 
69790 11V3 
SYVUOTIVG 
L¥OQOTIVG 
ooo liye 
G7 OONL iS 
YvVOO TIVA 
cvOOTIVS 
cCVUITIVG 
bzu0G tiva 
O7VCOTVS 
o€G01IV4 
SeO0ONIVves 
iCoOo Vd 


3d= 1d 

S71 eee 

(TO*V)/4= id 

d=) 

OCT OCI Udi Ni oae ee 
(J dJ2j 1395-1) 

OD Cae Os) 

Jd= Vd 

V/9=V9 

Oct *Sc1>Oul (yeaa 
(Jd/s) ives 2 

OCT *IECT Ge Ge vai 
OLE* S¥e* Gye HW Nao) ee 
J0+a>3 


(T°64¢°9 "8 st ° 14° 1 6a" 9 BS" l 6 ee ee ee 
N° SULS*® X3°V 0° So ACOs ee ees 


SolS* on 
(X2)SVWOIS=SULS 


C(cexedd-O°L)/(O0°TH(G1SLAS9ST* LI /(H/S) J eW2 GeO" <7 4) -—xe 
GO? *O00e *O0S (G=XaG Ei an 


T#X30LIT=As 0s 
NS=S 

OOc O1 saa 

O*Z/ (S=NS } +See 
O° O=u3 ICL 


O6¢*S6¢*SS5e (0°S2—-Ys1GL lee 


(S-NS) SU? =a 
(NOLS JAVWNOTS &H=NS 
Si” 31 Ge 
HexQtslA=Nns$ 


GOO"OTY “OL> (Age See 
(NE LS ud evo VaLNV AS Sits 


V/M=VM 

O6C*OLE* OTE (aaa) eese 
WZOx«(9-V)*«S*° OtWZ=EM 

(S xO °?}/ (00% Dx Td) =H! 70 


oat 
Oct 


Vicar 


Oye)! 


SCC 


G7e 


OSE 
OBE 


Ote2 


G62 


S6¢ 


cl? 
SOY 


OOTY 
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T8001IV4 
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